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Abstract 

We present a measurement of CP-violation parameters in the b — > sqq penguin transitions 
{q = s,u,d) based on a 78 fb~^ data sample collected at the T(4S') resonance with the Belle 
detector at the KEKB energy-asymmetric e+e~ collider. One neutral B meson is reconstructed 
in the (t)K^, K'^K~Kg, or rj'Kg decay channel, and the flavor of the accompanying B meson 
is identified from its decay products. CP violation parameters for each of the three modes are 
obtained from the asymmetries in the distributions of the proper-time intervals between the two 
B decays. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the standard model (SM), CP violation arises in weak interactions from an irreducible 
complex phase in the Kobayashi-Maskawa (KM) quark- mixing matrix []}]. In particular, the 
SM predicts CP-violating asymmetries in the time-dependent rates for and decays 
to a common CP eigenstate fcp 0- Recent measurements of the CP- violating parameter 
sin 201 by the Belle p, 0| and BaBar collaborations established CP violation in neutral 
B meson decays mediated by the b — >• ccs tree transition at a level consistent with KM 
expectations. 

Despite this success, many tests remain before one can conclude that the KM model 
provides a complete description. For example, the charmless decays B^ P° 
K~^K~Kg, and P° v'^s^ which are mediated by the b sss transition {B^ v'^s 
also receives b —>■ sdd and b suu penguin contributions) are potentially sensitive to 
new CP- violating phases from physics beyond the SM |7i]. SM contributions from the 
b — > VMS tree diagram are expected to be highly suppressed [1, 0| . Thus, the SM predicts 
that CP violation measurements in these charmless modes should yield sin20i to a good 
approximation. Consequently, a significant deviation in the time-dependent CP asymmetry 
in these modes from what is observed in 6 — > ccs decays would be evidence of a CP- violating 
phase not expected in the KM model. 

In the decay chain T(4S') B^B^ — > fcpfta.g, where one of the B mesons decays at 
time tcp to a final state fcp and the other decays at time ttag to a final state ft^c, that 
distinguishes between P° and P*^, the decay rate has a time dependence given by [10| 

ViAt) = — \l + q- Ssm{AmdAt) + Acos{AmdAt) }, (1) 

4:Tbo I l J J 

where tbo is the P° lifetime, Am^ is the mass difference between the two P° mass eigenstates. 
At = tcp — ttag, and the 6-flavor charge q = +1 (—1) when the tagging B meson is a P° 
{B^). The CP- violating parameters S and A are given by 

^_2Jm(A) |Ap-l 
\X\' + V ^=|AP + 1' 

where A is a complex parameter that depends on both the B^B^ mixing and on the am- 
phtudes for P° and P° decay to fcp- To a good approximation in the SM, |A| is equal to 
the absolute value of the ratio of the B^ ^ fcp to B^ — > fcp decay amplitudes. The SM 
predicts S = — ^/sin20i, where ^/ = 1) corresponds to CP-even (-odd) final states; 

and ^ = (or equivalently |A| = 1) for both b — > ccs and b — >• sss transitions. 

In this paper, we report the first measurement of CP asymmetries in P° (p^s 
K^K~K'g decays, and an improved measurement for P° v'^s decay based on a 
78 fb~^ data sample, which contains 85 million BB pairs. Data are collected with the Belle 
detector at the KEKB energy-asymmetric e~^e^ (3.5 on 8 GeV) collider |12| operating at 
the T(4S') resonance. At KEKB, the T(4S') is produced with a Lorentz boost of = 0.425 
nearly along the electron beamline (z). Since the B^ and B^ mesons are approximately at 
rest in the T(4S') center-of-mass system (cms). At can be determined from the displacement 
in z between the fcp and /tag decay vertices: At ~ {zcp — ^tag)//37c = Az/Pjc. 

The Belle detector jl^ is a large-solid-angle spectrometer that includes a three-layer sil- 
icon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), time-of-fiight (TOF) scintillation counters, and an 
electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect mesons and to identify muons (KLM). 
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We reconstruct decays to (pK^ and rj'K^ final states for = —1, and 5° K~^K~Kg 
decays tliat are a mixture of ^/ = +1 and —1. K^K~ pairs that are consistent with 
K^K~ decay are excluded from the K~^K~Kg sample. We find that the K~^K~Kg 

state is primarily = +1; the ^/ = +1 fraction is 1.04 ± 0.19(stat) ± 0.06(syst) In the 
following determination of S and we fix = +1. The intermediate meson states are 
reconstructed from the following decay chains: t]' ^ 7r"'"7r~)7 or t]' 7T~^TT^r]{^ 77), 

K^s 71+71-, and ^ K+R-. 

Candidate Kg —>■ tt+tt^ decays are oppositely charged track pairs that have an invariant 
mass within 15 (12) MeV/c^ of the nominal mass for the 5° (pK^ K+K~Kfj 
mode. The displacement of the tt+tt^ vertex from the nominal interaction point (IP) in the 
plane transverse to the positron beam axis (r-0 plane) is required to be greater than 0.1 cm 
and less than 20 cm. The direction of the combined pion-pair momentum in the r-0 plane 
is required to be within 0.2 radians of the direction defined by the IP and the displaced 
vertex. We select candidate (f) K+K~ decays requiring that the K+K~ invariant mass 
is within 10 MeV/c^ of the nominal meson mass, the meson momentum in the cms 
exceeds 2.0 GeV/c, and the K+K~ vertex is consistent with the IP. Since the meson selec- 
tion is effective in reducing background events, we impose only minimal kaon-identification 
requirements. For selection of non-resonant K+K~K^ candidates, more stringent kaon- 
identification requirements, which retain 86% of K+ at a 7% fake rate for vr^, are used. In 
this case, charged tracks that are positively identified as electrons or protons are excluded. 
In addition to the rejection of meson candidates, we reject K+K~ pairs that are consistent 
with K+K- or Xco ^ K+K' decay. 

For reconstructed B ^ fcp candidates, we identify B meson decays using the en- 
ergy difference /\E = E'^^ — -Ebeam beam-energy constrained mass Mbc = 
\f{E^^^^)^'^^Jp^^ , where E^^^-^ is the beam energy in the cms, and E]^^ and p^*^ are 
the cms energy and momentum of the reconstructed B candidate, respectively. The B me- 
son signal region is defined as 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and \ AE\ < 0.051 GeV for 
B^ 4>Kg or \ AE\ < 0.040 GeV for B^ —>■ K+K^Kg. In order to suppress background from 
the e+e~ QQ {q = u, d, s) and cc continuum, we form signal and background likelihood 
functions, £3 and £bg, from a set of variables that characterize the event topology 
We determine Cs from Monte Carlo (MC) and Cbg from data, and impose mode-dependent 
thresholds on the likelihood ratio Cs/i'^s + -^bg)- The numbers of reconstructed candidates 
are 59 and 230 for 5° ^ 0i^g and for 5° K+R-R^s, respectively. 

For B^ — > rj'Rg decay, we use the same selection criteria as those used in our previously 
published analysis |11|] if both of the charged pions in the rj' — »• 71+7t~ rj or the r]' — > p°7 decay 
have associated SVD hits. We also reconstruct events where only one of the charged pions 
has associated SVD hits. In this case, the requirement on the impact parameter is relaxed 
for the track without SVD hits, while a higher threshold is imposed on the likelihood ratio. 
The number of reconstructed B^ — > v'^^ candidates is 311. 

Charged leptons, kaons, pions, and A baryons that are not associated with the recon- 
structed fcp decay are used to identify the 6-fiavor of the accompanying B meson, which 
decays into /tag. Based on the measured properties of these tracks, two parameters, q and r, 
are assigned to each event. The first, g, has the discrete value -|-1 (—1) when the tag-side B 
meson is more likely to be a B^ (B^). The parameter r is an event-by-event MC-determined 
flavor-tagging dilution factor that ranges from r = for no flavor discrimination to r = 1 for 
an unambiguous flavor assignment. It is used only to sort data into six intervals of r, accord- 
ing to the estimated flavor purity. The wrong-tag probabilities for each of these intervals. 
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TABLE I: The numbers of reconstructed B fcp candidates used for S and A determination, 
Ncv, and the estimated signal purity in the Ai^-M^c signal region for each fcp mode. 



Mode 


if 




Purity 




-1 


53 




K+K-K^ 


+1(100%) 


191 






-1 


299 


0.49 ±0.05 



Wi {I = 1,6), which are used in the final fit, are determined directly from the data. Samples 
of decays to exclusively reconstructed self-tagging channels are utilized to obtain wi using 
time-dependent B^-B^ mixing: (A^^qf — No^y) / {Nq^ + A^^sf) = (1 — Stfz) cos(AmrfAt), where 
A^OF and A^sf are the numbers of opposite {B^W B^W) and same {B^W B°B^,WW) 
flavor events. The event fractions and wrong tag fractions for each r interval are described 
elsewhere 

The decay vertices of B^ mesons are reconstructed using tracks that have enough SVD 
hits: i.e. both z and r-0 hits in at least one SVD layer and at least one additional layer 
with a z hit, where the r-0 plane is perpendicular to the z axis. Each vertex position is 
required to be consistent with the IP proflle, which is determined run-by-run and smeared in 
the r-(j) plane by 21 fim to account for the B meson decay length. With these requirements, 
we are able to determine a vertex even with a single track. The vertex position for the fcp 
decay is reconstructed using charged kaons for B^ 0-^5 and B^ — >■ K~^K~Kg decays and 
using charged pions from or t]' decays for B^ — > f]'K^. The algorithm for the /tag vertex 
reconstruction is chosen to minimize the effect of long-lived particles, secondary vertices 
from charmed hadrons and a small fraction of poorly reconstructed tracks 3|. From all 



the charged tracks with associated SVD hits except those used for fcp, we select tracks 
with a position error in the z direction of less than 500 /im, and with an impact parameter 
with respect to the fcp vertex of less than 500 /im. Track pairs with opposite charges are 
removed if they form a Kg candidate with an invariant mass within ±15 MeV/c^ of the 
nominal mass. If the reduced associated with the /tag vertex exceeds 20, the track 
making the largest contribution is removed and the vertex is refltted. This procedure is 
repeated until an acceptable reduced is obtained. 

After flavor tagging and vertex reconstruction, we obtain the numbers of B^ ^ fcp can- 
didates, A^ev, listed in Table HI Figure [T] shows the Mbc distributions for the reconstructed 
B candidates that have AE values within the signal region. To assign an event-by-event 
signal probability for use in the maximum-likelihood flt of the CP-violating parameters, 
we determine event distribution functions in the AE-M]^c plane for both signal and back- 
ground. The signal distribution is modeled with a single two-dimensional Gaussian, where 
the widths are allowed to float in the flt to data. For the continuum background, we use 
a linear function for AE and the ARGUS parameterization for Mbc. We use events 
outside the signal region as well as a large MC sample to study the background compo- 
nents. The dominant background comes from continuum events. In addition, according 
to MC simulation, there is a non-negligible (~ 8%) contamination from BB background 
events in B^ v'^s iv' ~^ P^l)- The contributions from BB events are smaller for other 
decay modes. The contamination of K^K^Kg events in the (pKg sample (and vice versa) 
is also small and is treated as a source of systematic uncertainty. Finally, backgrounds from 
B^ fo{980)Kg decay, which has the opposite CP eigenvalue to 4>Kg, are found to be 
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FIG. 1: The beam-energy constrained mass distributions for — > (j)K^ (left), — > K~ Kg 
(center), and B^ v'^s (right) within the AE signal region. Sohd curves show the fit to signal plus 
background distributions, and dotted curves show the background contributions. The background 
for B^ — > rj'Kg decay includes an MC-estimated BB background component. 



negligible. 

We determine S and A for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for 
the signal distribution is given by Eq. ((T)) with q replaced by g(l — 2wi) to account for the 
effect of incorrect flavor assignment. The distribution is convolved with the proper-time 
interval resolution function Rsig{At), which takes into account the finite vertex resolution. 
It is formed by convolving four components: the detector resolutions for zcp and Ztag, the 
shift in the Ztag vertex position due to secondary tracks originating from charmed particle 



decays, and the kinematic approximation that the B mesons are at rest in the cms j3 . 
A small component of broad outliers in the Az distribution, caused by mis-reconstruction, 
is represented by a Gaussian function Poi(At). We determine twelve resolution parameters 
and the neutral- and charged-S lifetimes simultaneously from a fit to the At distributions 
of hadronic B decays and obtain an average At resolution of ~ 1.43 ps (rms). We determine 
the following likelihood value for each event: 



/oo r 
/sigPsig(At', q, wi)R,,^{Ati - At') 
-oo L 



+ (1 - /sig)Pbkg(At')^bkg(Ati - At') 



d(At') + /oiPoi(At, 



(3) 



where fo\ is the outlier fraction and fsig is the signal probability calculated as a function of 
AE and Mbc. Pbkg(At) is a PDF for background events, which dilutes the significance of 
CP violation in Eq. (^. It is modeled as a sum of exponential and prompt components, 
and is convolved with a sum of two Gaussians, -Rbkg, which represents a resolution function 
for the background. All parameters in Pbkg(At) and i?bkg are determined by the fit to the 
At distribution of a background-enhanced control sample i.e. events away from the 



AE-M^c signal region. We fix the r^o and Am^ at their world-average values jl7|. The 
only free parameters in the final fit are S and A, which are determined by maximizing the 
likelihood function 

L = l[P,{Atf,S,A) (4) 

i 

where the product is over all events. Table ITTl summarizes the results of the fit. The table 
shows the values of A and —C,fS, which, in the SM, is equal to sin20i. The first errors are 
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TABLE II: Results of the fits to the At distributions. The first errors are statistical and the second 
errors are systematic. The third error for the K~ mode arise from the uncertainty in the 
fraction of the CP-odd component. 



Mode 


-^fS (= sin2(/)i in the SM) 


^ (= in the SM) 




-0.73 ±0.64 ±0.22 


-0.56 ±0.41 ±0.16 


K+K-K^ 


±0.49 ± 0.43 ±0.1lt[5:i];j 


-0.40 ± 0.33 ±0. 10 j:[J:^|] 




±0.71 ± 0.37l[5:o6 


±0.26 ±0.22 ±0.03 




FIG. 2: The At distributions for (pK^ (left), 5° ^ K+K-R^ (center), and S° rj'K'^ 

(right) decays. The upper and lower plots are for q(^f = —1 and g^/ = ±1 candidates, respectively. 
The solid curves show the results of the global fits, and dashed curves show the background 
distributions. 



statistical and the second errors are systematic. The third error for the K^K~ mode 
arises from the uncertainty in the fraction of the CP-odd component Figure |21 shows 
the observed At distribution for g^j = — 1 (upper figure) and q^f = ±1 (lower figure) event 
samples for each fcp mode. Figure 121 shows the raw asymmetry in each At bin without 
background subtraction, which is defined by 

iV,e,.=_i±iV,5,.=+i' ^ ^ 

where iVg^^=+i(_i) is the number of observed candidates with qC,f = ±1(— 1). The curves 
show the results of the unbinned-maximum likelihood fit to the asymmetry distribution, 
— ^jiS sin(AmrfAt) — ^j^cos(AmrfAt). These are the first measurements of the CP violation 
parameters for — » dK^ and 5° K^K~K^ decays. The result for rj'K^ supersedes 
the previous result [ll|. We obtain values consistent with the present world average of 
sin 201 = ±0.734 ± 0.054 [l^ in 5° ^ K+K'K^ and r/'/sTg decays, while a 2.1a deviation 



is observed in (pKg decay. 

Fits to the same samples with the direct CP violation parameter A fixed at zero yield 
= -0.83±0.72(stat) for B^ (PK% -^fS = ±0.59 ± 0.47(stat) for 5° ^ K+R-R'^s, 
and —C,fS = ±0.77 ± 0.38(stat) for B^ v'^^s- -^^ ^ consistency check for the S term, we 
select the charged B meson decays B^ — > (pR^ and B^ —>■ rj'R^ and apply the same fit 
procedure. We obtain S = 0.05 ± 0.32(stat), A = 0.29 ± 0.24(stat) for B+ (pR+ decay 
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FIG. 3: The At asymmetry, A, in each bin for S° (jyK^ (left), 5° K+R-R'^ (center), and 
— > rj'Rg (right), respectively. The curves show the results of the unbinned-maximum likelihood 

fit. 



and S = -0.03 ± 0.20(stat), A = 0.05 ± 0.13(stat) for 5+ t]'K+ decay. Both results on 
the S term are consistent with no CP asymmetry, as expected. 

The largest source of systematic error for the 5° 0-^5 mode is the uncertainty in 
the signal fraction and the background At shape (±0.17 for S and ±0.14 for A in total) 
determined from the events in the sideband regions in the AE-M-^c plane. Other significant 
contributions come from uncertainties in the vertex reconstruction, the resolution function 
parameters, wrong tag fractions, r^o, and Am^. We add each contribution in quadrature 
to obtain the total systematic uncertainty. Systematic uncertainties from these sources are 
also examined for the other modes. We find that the largest uncertainties arise from the 
vertex reconstruction (±0.09 for S and ±0.08 for A) for the —>■ K^K~K^ mode, and 
from the resolution function parameters {^^^i for S) and the signal fraction (±0.02 for A) 
for the — > rj'K^ mode. 

In summary, we have performed the first measurement oi CP violation parameters in 
the B^ (p^s and K^K~Kg decays. We also provide an improved measurement for rj'Kg 
decay. These modes are dominated by the b sss transition and are sensitive to possible 
new CP-violating phases. Our results for P° — >■ rj'K^ and K^K~Kg are consistent with 
those obtained for P° — J/ipKg and other decays governed by the b ccs transition. A 
2.1(7 deviation is observed for P" — * 0-^5- 
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